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The quantum-mechanical A/'-body problem has deserved a tremendous amount of work since more than 50 years. 
Both numerical and analytical resolutions of that problem are highly nontrivial as soon as A/" > 2. Actually, even 
restricting to two-body systems, an analytical solution is generally not known although numerical computations are 
quite straightforward in this case. Because quantum A/'-body problems arise in all areas where quantum mechanics is 
involved, from atomic to hadronic physics, the amount of papers devoted to that topic is huge. We refer the reader 
to several textbooks in which many analytical results as well as useful references can be found [J-O ■ 
r^. Recently, we have proposed a new method, called auxiliary field method (AFM), to compute the bound states of 
1^ ' a given two-body Hamiltonian 0-01 • This method is based on auxiliary (einbein) fields and can lead to analytical 
, approximate closed formulas for various central interactions, with both nonrelativistic or relativistic kinematics. The 
dependence of mass formulas on parameters of the Hamiltonian can generally be determined from a simple dimensional 
analysis. But one of the interesting features of the AFM is to predict, at least approximately, the dependence on 
quantum numbers. 

' The purpose of the present paper is to show that the AFM can be successfully applied to find approximate analytical 
mass formulas for general A/'-body Hamiltonians of the form 



^ = E \/p' + + E ^'(1^^ - ^1) + E i^^.(in-^.i)' (1) 
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with X^i^iP^ = P = 0, where and Pi are respectively the position and momentum of particle i with a mass m^, 
I and where i? is a global variable to be specified later. In the following, we will use the natural units h = c = I. 
•• i We consider a relativistic kinetic energy and allow for both pairwise and one-body interactions. Such a Hamiltonian 
^ ' being not Lorentz invariant, it is often called "semirelativistic" in the literature, see e.g. [10, 11]. Particular systems 

. exhibiting one-body potentials of the form appearing in Eq. ([T]) will be studied below. 
}^ ' Our paper is organized as follows. We begin by showing how the ground-state energy of a A/'-body Hamiltonian with 
5^ 1 two-body interactions can be estimated from the corresponding two-body Hamiltonian in Sec. |TT1 Then, we focus on 
the A/'-body harmonic oscillator in Sec. IHII This problem is particularly interesting since it can be fully analytically 
solved without resorting to the AFM in the nonrelativistic case, while the AFM is needed in the semirelativistic case. 
The harmonic oscillator is the starting point of the general application of the AFM to the generic Hamiltonian ([T]) . The 
procedure is discussed in Sec. IIVI followed by several examples: power-law potentials in Sec. |Vl baryonic and atomic 
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systems in Sees. I VII and I VIII The case for which only two-body interactions are present is discussed in Sec. IVIIII 
Finally, an application of our formulas is given in Sec. IIXI and conclusions are drawn in Sec. |Xl 

II. TWO-BODY REDUCTION 

Finding exact results for two-body Hamiltonians is obviously a less problematic task than in the A/'-body case, 
not only analytically but also numerically. That is why it is of interest to try to reduce A/'-body Hamiltonians to 
equivalent two-body problems. That idea is already present in early works about nuclear physics such as Ref. p^. 
Further references about that technique can be found in Ref. [isj. In the same spirit, we show in this section how 
the AFM can be applied to express the ground-state energy of a A/'-body system in terms of the one of a two-body 
system. 

A. General results 

Several analytical approximate energy formulas have already been found by applying the AFM to one- or two- 
body Hamiltonians, see Refs. 043- These formulas can actually be used to compute the ground-state mass of a 
system of M identical particles interacting via pairwise interactions. Let us start from the A/'-body spinless Salpeter 
Hamiltonian 

U A/" 

= J^y^pf + m2+ ^ V{n,), where = |n - r,|, (2) 

i=l i<j = l 

which is a particular case of Hamiltonian ([T|). We denote M'^-'^^m) the exact ground state mass of the system governed 
by the Hamiltonian ^ ; we want an approximate analytical expression of this quantity. 

Finding analytical solutions for a wide range of potentials would be hopeless because of the square roots of the form 
•y/p^ + irfi appearing in the kinetic term. Nevertheless, as already pointed out in Ref. the operator inequality 

V^^<^^^f^, (3) 

where // is a real number, can be used to obtain upper bounds on the exact semirelativistic Hamiltonian. Analytical 
results can then be hoped since the introduction of the variational parameter jj, leads to an apparently nonrelativistic 
kinetic energy. Such a technique has been further developed in Refs. [13, [H, [ill, where explicit example are also 
given. In our framework, the auxiliary fields, denoted as Vj^ play the role of the parameter /i. When introduced in 
Hamiltonian ([2]) in order to get rid of the square roots + m? , one obtains 



We refer the reader to Ref. [1] for a detailed discussion about this procedure, but recall here the main ideas 
for completeness. The auxiliary fields are a priori operators that can be eliminated by solving the constraints 
= 0, which are the Euler-Lagrange equation for the field Vi. It can then be checked that 

Vi — M \J pI -f m2/2 and that H'^^^(vj) = H^-^\ Both Hamiltonians ([2]) and (Hj) are equivalent up to a proper 
elimination of the auxiliary fields. The computations are however considerably simplified if one makes the approxi- 
mation that the auxiliary fields are real numbers. Then, Hamiltonian ^ becomes formally a nonrelativistic one, and 
the auxiliary fields have to be seen as variational parameters, whose optimal values, Vj.o^ are eventually obtained by 
minimizing the eigenvalues of (HJ, denoted as M^-'^^vjim). Provided that Hamiltonian (|4]) is analytically solvable or 
that an upper bound of its eigenvalues can be found, the inequality ([3]) says that M^^^iyj^; m) is an upper bound of 
the exact mass M(^)(m), that IS 

M^-^\m) < M(-^)(i/j-o;to) < M^-^^iiyf,m). (5) 

The optimal values Vj^ depend on quantum numbers of the state considered. If a lower bound of the eigenvalues 
of (01) is found however, no conclusion can be drawn about the relative position of this last bound and M^-'^\m). 
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Prom the AFM, we know that the optimal vahies vj.o should be close to {Dj) The ground state being the most 

symmetrical state, it is reasonable to assume that Vj = v hy symmetry, i.e. that the total kinetic energy is equally 
distributed in average between the different particles. This will shown in Sec. IIV A2I With such a choice, the AFM 
can be applied with a single auxiliary field v instead of M auxiliary fields Vj , a much simpler situation. Consequently 
one is led to 



/■2^2 «r ^f 



i—1 i<j—l 

Let us denote by M'^^\v; m) the ground state mass of the simpler Hamiltonian ^ and i/q the value that minimizes 
this quantity. Due to the fact that the mass depends now on a single auxiliary field, the expression of M'^-'^\i^o;m) 
is much easily obtained than M(-^^(z^j^o; w)- As M^-'^^i^; m) = M^^)(z^j — v]m), we can write 

M(^)(m) < M^^^\vo-m) < M^^^\v;m). (7) 

Let us consider the two-body nonrelativistic problem of two particles of mass v interacting through the potential 
V{r). The corresponding Hamiltonian is thus /i^ + V{r) whose ground state energy is e(i^). It is convenient to use 
the following identity Q 



^^Hp^-^ E 4+ EpO 



(8) 



where 77^ — ^' is the relative momentum of the pair (ij). Since we are working in a frame where Pi = Oj 

the second term in the right-hand side of ([8]) cancels and Hamiltonian ([4]) becomes exactly 



i<j=l L 



V 



(9) 



Let us assume that we know an analytical expression of t{v). Since we are dealing with a system of M particles in 
their ground state, and since the Hamiltonian (fSl) is separable, it can be shown that a lower bound on the ground-state 
mass M^^-' (v) of Hamiltonian ([9]) is given by Q 

M(^)(.;m) > M(^)(.;m) = ^ ^ , ^^^^ ^ "^Kj^,). (10) 

Now, let us consider the semirelativistic version of the two-body problem, whose Hamiltonian is given by i?*-^) — 
a \Jp'^ -\- rn? -\- Vir). In principle the two-body case implies cr = 2. However, for further convenience, a is chosen as 
an arbitrary positive real number. The exact ground state mass of this Hamiltonian is denoted M'^^)(m, cr). 

Introducing the auxiliary field /z, one can use AFM to have an approximate expression of this quantity. Explicitly 
the auxiliary Hamiltonian is 

9 2 2 

(^) ^ _ + ^ + ^ + v{t). (11) 
The corresponding ground state Af '^^^ {p; m, cr) is given by 

M(2)(^;m,a) = — + ^-Ke(A.). (12) 

One has the following inequalities 

M(2)(m,o-) < A/(2)(^o;TO,cr) < M'-^\fj.]m,a), (13) 
where fj,o is the value minimizing Af '^'(/i; m, cr). Comparing ([TU|) to (|12p . one can easily show that 
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Consequently, we arrive at the following series of inequalities M^^\vo;m) > M^^\vo]m) ^ {N{N - 
l)/2)M(2)(i/o;m',a') > (AA(A/' - 1)/2)M(2)(^o; m', a') > {U{U - I) /2)M^^\m' ,o'). 
The previous inequalities and inequalities (O are summarized by 

M^^\m) < M^^\vn;m), (15) 
^^M_AM^^\m',a') < M^^\iyo;m). (16) 

If M^-^^ {vq; to) represents a good approximation of the A/'-body energy, it seems reasonable to assume that the values 
of the left hand parts should be rather close. We arrive at the final approximate formula for the A/'-body ground-state 
mass in terms of the two-body ground-sate mass 

This means that all the ground state mass formulas that can be obtained in the two-body case can be generalized to 
the ground state of a A/'-particle system thanks to p7)) . It is worth pointing out that the result (fT7)) has been firstly 
conjectured in Ref. 17[ and explicitly proved for A/" 3, 4 in this last reference. It is an interesting check of the 



AFM that we recover this previously established result. For completeness, we mention that other lower bounds can 
be found for example in Refs. [isl fioj. 

Let us emphasize an important point. In this section the auxiliary field v was introduced to deal only with the 
worrying square root appearing in the kinetic energy operator. This means that, in order to apply (I17p . one must be 
able to obtain the exact eigenvalues of the Hamiltonian H'^^^v) [see Eq. pTjl ] with the given potential V{r). In that 
case, the spectrum obtained is an upper bound of the exact one. Indeed, very rare are the situations for which this 
procedure can be achieved. More frequently, one has only approximate expressions for e(y)\ in this respect the AFM 
can be invoked once more, but at the price of introducing supplementary auxiliary fields [J|. Then, an upper bound 
of the spectrum is obtained if an upper bound of e{v) can be found, and no conclusion can be drawn if a lower bound 
of e{v) is obtained. This approach, dealing with the kinetic energy and potential operators on equal footing, will be 
discussed in further sections. 



B. An example 

As an example of the above considerations, let us consider the case of A/" identical particles interacting via a Coulomb 
potential — . Applying the AFM, the ground-state mass of Hamiltonian o^p^ -\- nn? — b/r is given by Q 



Mg,,{m,a) = <7m^ll~—^ . (18) 



In this last relation, Qc ^ n -\- 1 + \ a priori^ but it can also be replaced by a more complicated function of the form 
Qc = a(h) n + l + (3{b), with a and /3 continuous functions of b, giving better results Q. Equation (|T7)) leads then to 

With a relativistic kinematics, the strength of the Coulomb potential must not be too strong otherwise the Hamiltonian 
is not bounded from below and the system collapses. It is easy to notice that the critical value, b^f, for the AA-body 
problem is related to the critical value, 62, for the two-body problem through 



Thus it decreases as the number of interacting pairs increases. We also recover a result given in Ref. |20|,|21j: At large 
A/", the stability of the system is guaranteed if 5 = b/M, where 6 is a constant below some critical value. 
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III. THE A^-BODY HARMONIC OSCILLATOR 

The case of a quadratic potential, i.e. the A/'-body harmonic oscillator, is practically the only one for which an 
exact solution is reachable. For that reason, it is the basic ingredient of our method when dealing with AA-body 
Hamiltonians. The A/'-body harmonic oscillator problem has been widely studied in the literature, especially in the 
nonrelativistic case. A first analytical solution of that problem can be found in Ref. (22j together with an application 
to nuclear physics, while the exact spectrum of a nonrelativistic A/'-body harmonic oscillator with pairwise interactions 
has been given in Ref. 23]. The interested reader will find a detailed list of references concerning that topic in Sec. Ill 



of Ref. [2JI, as well as a discussion of the changes caused by the consideration of fermions rather than bosons. In the 
semirelativistic case, upper and lower bounds are found in Ref. [iJl for a A/'-body harmonic oscillator with pairwise 
interactions. 

Before making explicit computations, it is worth mentioning the new results reported in this Section. First, we 
extend previous results in the nonrelativistic case by showing that exact eigenenergies can be found when both one- 
and two-body harmonic interactions are present [see Hamiltonian (1261) ]. Second, we give an explicit mass formula in 
the semirelativistc case with one- and two-body interactions, which is an upper bound of the exact mass formula. To 
our knowledge, that formula has never been obtained before. 

A. Jacobi coordinates 

Usually the many-body problem is treated starting with the shell-model or variants. In this method, the degrees of 
freedom are simply the positions of the various particles (i = 1, . . . , A/'). But, in this case, the motion of the center 
of mass is not separated correctly and this leads to spurious components which spoil the results. On the contrary, in 
the few-body problem, one introduces new degrees of freedom which allow to correct exactly this important drawback. 
The price to pay is that the resulting method becomes more and more difficult to be applied when the number of 
particles increases. In practice, one is limited to values less than 10. Fortunately, this method can be used in the case 
of the harmonic oscillator, whatever the value of Af, as we will see below. 

Among the various possible new degrees of freedom, the Jacobi coordinates are of common use. First, we choose a 
reference mass m (it can be that of a particle or the total mass of the system for example) and define the dimensionless 
quantities 

i 

ai=TOz/m, Q;i2...i = ^afc and a^ai2...Af- (21) 

fc=i 

Then, the standard Jacobi coordinates, be expressed as 

Efe=iafcrfc ^ i = l,...,JV -1, and xj^ = R = - J" akVk. (22) 



Xi = — — n+1 

ai2...j a 



k=^^ 



For convenience, the global coordinate R is relegated as the last Jacobi coordinate xj^f whereas Xi {i — I, . . . , Af — I) 
represents the vector joining the nonrelativistic center of mass of the first i particles to the particle i -I- 1. In the 
case of nonrelativistic kinematics or in the case of a system composed of identical particles, R is the center of mass 
coordinate. By using the matrix notation Xi = ^v'^'j^ O'^^ ^^'^ to the following definition of the [/-matrix 

U^j = if j <i (23a) 
ai2...i 

Uu+i = -1 (23b) 

U^J =0 if j > i + (23c) 

It is not difficult to calculate the inverse matrix B = U^^ whose matrix elements are given by 

Bki = -^^i±i_ if k<l<N (24a) 
OL12...1+1 

Bi+ii = -^1^ if KM (24b) 
ai2...i+i 

Bki ^ if fc > / -I- 1 (24c) 

BkAT = 1 Vfc. (24d) 
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Denoting pi and tt^ the conjugate momenta associated respectively to and Xi, it is easy to prove that 

A/" Af 

TTi = ^ BjiPj and = ^ U■i^^Zj. (25) 
TTj\f — P is the total momentum of the system. In any case, we work in a frame for which TVj\f — 0. 

B. Nonrelativistic general case 

Let us start with the most general harmonic-oscillator-like Hamiltonian, corresponding to Af particles of arbitrary 
masses, with one-body and two-body quadratic potentials with arbitrary positive string constants. The Hamiltonian 
looks like 

= E ^ + E - + E -"-of- (26) 

i—1 i—1 i<j — l 

The kinetic energy operator, expressed in terms of Jacobi variables, allows the correct separation of the center of 
mass motion and appears decoupled in the various variables (this is in fact the justification of the form of Jacobi 
coordinates) 

T = y^ = —+y ^nf, (27) 

i—1 i—1 

where rrit = am = mi + m2 + ■ ■ ■ mj\f is the total mass of the system. The kinematical quantity Ai is calculated as 

/ \ 1/2 

A, = ( ""■■■'^^ ) . (28) 
\ai+iai2...ij 

For further convenience, it is judicious to switch from the standard Jacobi variables to conjugate renormalized Jacobi 
coordinates defined as (i = 1, . . . , A/" — 1) 

yi = Y'' Pi = KT^f (29) 

Working in the center of mass frame {P — 0) and using these new variables, the kinetic energy operator has a very 
simple form 



T-^'t.l (30) 

With the condition P = 0, it is clear that the case M — 1 cannot be treated by the following formulas. Results for 
this case are given in Ref. 

With these new variables the one-body operator is written 

Vi^Y^ k,{n -Rf= Pi"^yi ■ (31) 

i—1 /,m— 1 

where the symmetric positive-definite matrix F is defined by 

AT 

Flm — A;Am ^ hBuBim, (32) 
1=1 

the matrix B being given by (|24p. In the very same way, the two-body operator is written 

JV A/"-! 

= ^ kij{ri - TjY = ^ GhnyrUm, (33) 

i<j—l l,m—l 
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where the symmetric positive-definite matrix G is defined by 



Glm = ^l^m, ^ kij{Bii — Bji){Birn — Bjm)- (34) 
i<j=l 



Introducing the matrix J = F + G, the total potential y = Vi + V2 is expressed as 



V = ^ JlmVl ■ Vm- (35) 
l.m—1 

For arbitrary masses or/and string constants there is no reason why the matrix J should be diagonal. However, this 
matrix being a symmetric positive-definite matrix, it can be diagonalized with help of a unitary matrix (in fact an 
orthogonal one since all the quantities are real). Thus (O denotes the transposed matrix of O) 

J = Q-^DO, with Q-^ = 6. (36) 

The elements of the diagonal matrix D are all positive and are chosen under the form di = miJf/i. 
The last step is an ultimate change of conjugate variables 

A/"-! N-l 

zi^^ OijVj and cr; = ^ Oijpj. (37) 

Expressed with these new variables, the original Hamiltonian (j26l) appears to be the sum of A/"— 1 decoupled harmonic 
oscillators 



(T ; 1 09 

— ^ -I- -mujfzf 
2m 2 * * 



(38) 



4=1 

and consequently the energy of the system is given by 

Sho - -^k + 3/2), (39) 

i=l 

where Ui and U are respectively the radial and orbital quantum numbers associated to the coordinate Zi . The problem 
is now completely solved. Moreover, this result is valid for any excited state. Even if the expression of uji is in general 
not analytical for A/" > 5, the result (15^1) is exact and could be calculated with a high accuracy. Notice that formula ([M)) 
extends previous results (23l. [25j. where an equivalent mass formula is obtained in the case fc^ = 0. 



C. The 3-body harmonic oscillator 

The application of the general formulas to the special case A/" = 3 is very interesting: A very accurate numerical 
solution of the 3-body problem is now commonly obtainable and the harmonic oscillator is an obvious test of the 
codes. Our formulas could then be compared to numerical solutions relative to the 3-body problem. In this case the 
J-matrix is a 2 by 2 matrix so that the eigenenergies can easily be obtained. Moreover the original Hamiltonian does 
not favor a given particle and the eigenenergies (j39p must reflect this property; in contrast, all the formulas given in 
the previous section rely on the definition of the Jacobi coordinates which give a particular status to particle 3 with 
respect to particles 1 and 2. It is then very instructive to check this symmetry property in the analytical expression 
of the energy in this particular situation. 

Applying the method given previously, the result is obtained after a long but straightforward calculation. As before 
we introduce the reference mass m, the relative masses = nii/m, the dimensionless quantities — ai + aj, 
a = ai +02-1-0:3, and the dynamical quantities kij = ki + kj (not to be confused with kij). Defining the symmetric 
quantities 

s — kia2a3a23 + fc2aia3ai3 + ^30102012 + a(fci2a3ai2 + fci3Q!2ai3 + fc23aia23), (40a) 
r = kik2al + kiksal + k2k3al + a'^ {ki2ki3 + ki3k23 + ki2k23) + 

ki2iki2al + ^3012) + ^13(^13^2 + ^2^13) + ^23(^23"? + kiaj^), (40b) 
6 = a/s^ — 4aiQ;2Q!3Q;r, (40c) 
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the eigenenergies of the system read 



ho 



1 


1/2 







(s + 5)i/2(2ni + h+ 3/2) + (s - <5)i/2(2n2 + ^2 + 3/2)1 . (41) 



They effectively present the property of not favoring one of the particles. Although this is not explicit, they do not 
depend on the reference mass m anymore. 

D. Case of identical particles 

It is of interest to rewrite the solution in the case where all the particles are identical, implying that they have the 
same mass rrii — m and the same string constants ki = fc, kij = k. Hamiltonian (|26p indeed reads in this case 

^1^° = E ^ + fcE(^'^ ^ + ~^ E - ^^)'- (42) 

i—l i—1 i<j—l 

Then, it is easy to see that the J matrix is already diagonal from the very beginning so that its eigenvalues are 
analytically known (they all read k + Mk) and, consequently, the eigenenergies of the system are also analytical. 
Explicitly they are given by {P — 0) 



^ho = \/-(fc+AAfc) g, (43) 

" TO 



where Q is the total principal number 



Q^J2{2n, + h) + -{N--l). (44) 

i=l 

One can check that Eq. (|4T|) reduces to Eq. (|43)) with J\f — 3 for identical particles. 

Wave function with good global quantum numbers and various symmetries can be built by appropriate linear 
combinations of states characterized by the same value of the global quantum number Q [26j | . All symmetrized states 
have then the same energy than the nonsymmetrized states. For spatial wave functions completely symmetrical, the 
ground state is obtained for the values rii = k = Wi, so that the principal quantum number is simply 

QsGS = ^(AA-l). (45) 

For mixed symmetry or completely antisymmetrical spatial wave functions, the situation is much more involved. An 
estimation of the ground state energy for the completely antisymmetrical case can be computed by choosing different 
values for the quantum numbers and piling the states (d identical values of the same Uj and Ij per state in order 
to take care of a possible spin-isospin-color degeneracy) up to the Fermi level. By considering only particle number 
insuring a saturated Fermi level (closed shell), one obtains 

QAGS = ^(AA-l)(i?/ + 2), (46) 
where Bf = 2nf + lf, the band number of the Fermi level, is the real positive solution of 



Asymptotically, we have 



N-l = ^{Bf + l){Bj + 2){B^+3). (47) 



81 \ 1/3^^4/3 
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E. Semirelativistic case 



Let us now consider the semirelativistic generalization of (|42p , namely 



A/" 



AT M 

i=l i<cj—l 



(49) 



Auxiliary fields, /x^, can be introduced in a similar way as it has been done in the previous section. One obtains 



^ r 2,2 
2 



1=1 



2^l, 



i—1 — l 



(50) 



Provided the total wave function is completely (anti)symmetrical, the optimal real auxiliary fields are all equals. A 
proof is given in the following section. Setting /i^ — /i, we obtain 



AT 



Af 



2fi 



(51) 



i<j=i 



whose mass formula is directly obtained from (j43p provided P = 

^2 ^ 



A'/ho(M) = ^ ( M + — ) + J-{k+JVk) Q 



(52) 



The remaining auxiliary field has to be eliminated by solving 9pAfho(/i)| ^ =0. Defining 



2m 2 - -s^ 

: X'^ and Y = — — 



2N' 



■2 \ 2/3 



(fc + A/'fc)g2 



the minimization condition is given by 



The only positive root of this last equation is given by 



X = G^{Y) = l^^) + yA/,My)-v{Y), 



with 



v{Y) = (^2+ V4 + y3 
The mass formula ((52)) finally becomes 

3 mt 



1/3 



Mho(Ato) 



Y2G-{Yy 



(4G_(r) + y) = 



-y (2 + V4 + r3 



2mt 



-1/3 



/3Y 



G-iY) 



+ G_(y)2 



(53) 



(54) 



(55) 



(56) 



(57) 



Using the results of Ref. [9|, it can be shown that these masses are upper bounds of the eigenvalues of H^o- Notice 
the simple ultrarelativistic limit 



lim Mho(/io) = t: [2Mik+Mk)Q 



211/3 



(58) 



It can be checked that formulas ([57)1 and ([55)1 reduce to the two-body relations obtained in Ref. ior N' — 2 and 
k = 0. We point out that the notation G_ {Y) has been chosen in order to keep the same symbols as in our previous 
works [EQ ■ 
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IV. AUXILIARY FIELD METHOD 
A. Main mass formula 

1. General case 



The AFM, introduced in Refs. 0, i], can be straightforwardly generalized to the case of A/'-body Hamiltonians 
of the form ([1]). The basic idea is to introduce auxiliary fields so that this Hamiltonian is formally replaced by a 
Hamiltonian for which an analytical solution can eventually be found. The auxiliary fields are denoted as fii, Vi, and 



Vij, and are introduced as follows 



i=l 



2 2fii 2fii 

+ [^'J + ^^dhi'^^j)) - % PiI^A'^^J))] . (59) 

i<j=l 

where (notice the definition of ^ \ri\) 

ri = In - R\, rij = |r,; - rj\, (60) 

and where 



where the prime denotes the derivative with respect to the argument. It is sufficient to set fik = nik Vfc to treat 
a nonrelativistic kinematics. We stress that, in order for the AFM to apply, it is assumed that the mass spectrum 
of ([55]) is analytically computable. In practice, this is only possible for the choice P{x) — P{x) = under some 
conditions as explained in Sec. IIIII 

The reason of such a definition is that Hamiltonians ([1]) and (|59p are equivalent if the auxiliary fields are properly 
eliminated. It can indeed be checked that 



5^j,^H{iik,Vk,Vki)\^^^f,^=^ => fM = ^pf + m1, (62) 
Siy,H{^ik,Vk,Vki)\^,.^0^ = Vi^Kiin), (63) 
^iyijH{^ik,Vk,Vki)\^^.^l,.^ ^ v^j ^ K^j{rij), (64) 

and finally that H^fik, Vk, vu) — H- 

The approximation underlying the AFM is now that, as we have done in previous sections, the auxiliary fields will 
be seen as real variational parameters. Let us note 

H{iikfl,Vkfl,^kifi)Wo) =^o|'/'o), (65) 

where /^fe.o, Vk.o, Vkifl are the real values of the auxiliary fields such that Mq is extremal. \lpq) is an eigenstate with 
fixed quantum numbers and a given symmetry. Using the same mathematical techniques than in Refs. 
can show that 



^ifc,o = \l {Vi)\p'i + m^tWa) with P = 0, (66) 

{(po\P{ri)\LpQ) = P{rifi) with r^^o = /^(i^i^o), (67) 

{ipo\P{rij)\LpQ) = P{rijfl) with nj^o = hjii^ij.o)- (68) 

It is also expected that ((^olAil'^o) ~ /^i,Oj {^o\i'i\^o) ~ ^^i.o, and {ipQ\i>ij\ipo) w i^ij,Q- Let us rewrite the Hamilto- 
nian (|59p under the form 

H{nk,Vk,Vki) ^T{^ik) + 'YUi{iyi,n) + ^ Uij{Pi,ri.j), (69) 

2=1 i<j = l 
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where T{^k) stands for the kinetic part. One can also check that 

Udi^uo, nfl) = T/(r,,o), {/;(m,,o, r,,o) = V'{n^o), (70) 
Uij{Pijfi,rij^o) = V{rijfi), U'ij{yijfi,rij^Q) = V'{rijfi). (71) 

This means that the approximate potential Ui {Uij) and the corresponding genuine potential Vi {Vij) are tangent 
at, at least, one point. If the following conditions Uiiyifi^ri) > V{ri) and Uij{vij^Q,rij) > V{rij) arc fulfilled Vi and 
Vj and for all values of the radial argument, Mq is an upper bound of an eigenstate oi H If the kinematics is 
nonrelativistic {T = T{^k = "m-k)) and if the following conditions Ui{vifi,ri) < V{ri) and Uij{vijfi,rij) < V{rij) are 
fulfilled Vi and \/j and for all values of the radial argument, Mq is a lower bound of an eigenstate of H. This restriction 
about T comes from the fact that the replacement of the genuine relativistic kinetic operator by the form T(/ife) yields 
upper bounds of the eigenvalues of the genuine Hamiltonian. In the other cases, it is not possible to obtain a relevant 
information about the position of Mq [3, |^ . 

At this stage, approximate numerical solutions of the A/'-body problem can be easily computed. First, the choice 
P{x) = P{x) = allows a precise determination of the eigenvalues of Hamiltonian ([5^ [see Sec. Ull]. Second, the 
extremization of the eigenvalues with respect to the real auxiliary fields is a classical numerical problem which can be 
solved with a high accuracy. 



2. Identical particles 



To obtain analytical closed solutions of the eigenvalue and extremization problems associated to Hamiltonian (j59p , 
it is necessary to simplify the system. First of all, we will only consider systems with identical particles, that is with 
rrii = m. In this case, it is reasonable to consider identical interactions between them, namely Vk{x) — V{x) and 
Vki{x) = V{x). This means that Kk{x) = K{x), Ik{x) = I{x), Kki{x) = K{x), and Iki{x) = I{x). Let us denote Pij 
the permutation operator exchanging particles i and j. We can write Pij|(^o) = ilv'o) if this state \lpq) is completely 
(anti)symmetrical. Then 



P{rjfi) = (¥'o|-P(r,)|(^o) - {^a\h,P{n)h,\^o) = {^o\P{n)\^o) = P(r,,o)- 



(72) 



If P{x) is monotonic, which is always the case in practice, then Vj^ — ri^. Finally, if I{x) is invertible, which must 
be the case to solve the problem, we have Vi^o = Klri^g) = K{rj^Q) — Vj.o- So all optimal values are the same. 
Using the same reasoning, we can draw the same conclusion for other auxiliary fields. 

Under these conditions, we can set = fi, Vi = v, and = P in the expression of the eigenenergies. Relations (j62p - 
(IMl) clearly show that they are related to the kinetic (/ij and potential energies (fi, Dij) in the system. We are thus 
led to replace the original Hamiltonian (1591) by the the following simpler Hamiltonian, which now depends on only 3 
auxiliary fields 



m 



-N[V{I{v))-vI{vf 



U{U - 1) 



[v[m)-Di{Df 



M 1 14 j\r 

i=l i=l 



(73) 



i<j=i 



which, by virtue of ([5^ . has the following mass spectrum {P = 0) 



+ 



m 



-M[V{I{u)) - vl{vf 



N{Af - 1) 



[V{i{9))-9i{uf 



(74) 



which is characterized by a high degeneracy due to the form of Q. The eigenstates of i?(/i, v, v) are built with harmonic 
oscillator states. As mentioned above, states with the same energy can be coupled to good quantum numbers and a 
given symmetry. 

The last step needed to get the final mass formula is to find the optimal values /^Oj and Pq from the extremization 
conditions 



9^Af(Ai,z^,i^)| =0, d^M{n,v,v)\ =Q, di,M{n,v,v)\-__ 



= 0. 



(75) 
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Writing explicitly Eqs. (I75|) and after some algebra, we arrive at the mass formula 



where the auxiliary fields are solutions of 



Mo 



m 
Mo 



MoAA2 



1/2 



1/4 



{U -lY^JiQ{vo+^^vo) 



1/4 



(76) 

(77) 
(78) 
(79) 



At this stage, functions / and / are not known. They depend on the specific forms of V and V . Optimal values of 
/io, va, and are no longer known. They also depend on the state considered through the variable Q. Nevertheless, 
formula (|76p makes clearly appear the mean-field nature of the AFM, because /io can be interpreted as the average 
kinetic energy of one particle, while V{I{vq)) and V{I{v{))) can be respectively seen as the average potential energy 
of one particle in the potential V{r) and of a pair in the potential V{r). 



3. General remarks 



A general remark about the AFM can now be done: The analytical results obtained for a given potential V{r) 
approximated by a potential P{r) can be used as a starting point for finding analytical results for an other potential 
VF(r), approximated this time by V{r). This considerably enlarges the domain of applicability of this method. 

Notice that the idea of rewriting a AA-body Hamiltonian wi th p airwise interactions of the form g{rfj) as a Af- 
body harmonic oscillator has already been investigated in Ref. [ll| within the framework of envelope theory. That 
method shares many similarities with the AFM, as shown in Ref. [7|, but the results presented hereafter have, to our 
knowledge, never been obtained so far. An interesting result coming from the comparison between AFM and envelope 
theory is that it allows to gain informations about the (anti)variational nature of the mass formulas obtained. Let 
us indeed assume a A/'-body Hamiltonian with pairwise interactions of the form Vij{rij) — g{rfj), where g{x) has a 
definite convexity. Then, provided that the kinetic term is exactly taken into account (i.e. in the nonrelativistic case), 
the energy spectrum obtained is an upper (lower) bound of the exact spectrum if g" < (> 0). For semirelativistic 
Hamiltonians however, an upper bound of the kinetic term is obtained because of the introduction of auxiliary fields. 
Consequently, in the semirelativistic case, the AFM leads thus to upper bounds of the exact spectrum if g" < 0. If 
q" > no conclusion can be drawn. 



B. Simplified formulations 



The four equations ([7S|) - (|7^ give the AFM approximation for the mass spectrum of a A/'-body system where the 
constituent particles are identical. They can be even simplified by defining 

Xo = \/2/io(i^o+AAPo). (80) 

We have indeed 

where the remaining auxiliary field is a solution of 

^0 - ^\l^' + JfXo 

and where K and K are defined by ((6T|) . This condition ensures that M{Xo) is extremal. The main technical 
difficulty is to find a closed analytical form for M{Xq). Fortunately, as we shall see in the next sections, this task can 




2Q 



(AA-l)Xo 



(82) 
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be achieved in several cases of interest. It is worth mentioning that the AFM results (I5T|) and (|82p are valid not only 
for the ground state but also for excited states. This is a serious advantage with respect to the two-body reduction 
presented in Sec. [TTl Identifying Eqs. (|76|) and (|8T|) . it appears that 



A.o-\/m2 + ^. (83) 



Moreover, the quantities /(t'o) and I{i'o) are directly obtained as functions of Xq. 



Ultrarelativistic and nonrelativistic limits 



For particles with null mass (ultrarelativistic limit), one obtains even simpler formulas by simply setting m = in 
formulas ^ and ([5^ . 

In the nonrelativistic limit, the auxiliary field fi tends towards m. In this case also the various formulas look simpler 
but cautions must be taken in the limit. Explicitly, Eq. (|8ip reduces to 



where the remaining auxiliary field is a solution of 



X^ = 2m 



AfXo \y{Af-l)Xa 



(85) 



V. POWER-LAW POTENTIALS 
A. General results 

The first explicit example that will be considered below is the case of power-law potentials, i.e. the Hamiltonian 

A/" ^^ A/" 

i/pi = ^ypf+^+asgn(A)^r,^+5sgn(r;) ^ 4, (86) 

2—1 2—1 i<j—l 

where sgn(p 7^ 0) = \p\/ p and where A, 77 > —1 (in the nonrelativistic case, one can consider A, > —2). When only 
a one-body or a two-body interaction is present, parameters a or 6 must be positive. If both types of potentials are 
present, there are less constraints on the sign provided that a bound state could exist. Following the definitions (j61L 
it is readily computed that 



Then, by defining 



Kix) = ^x'-', k{x)^'^x^~\ (87) 



2-A 2-77 



2Q 



Eqs. (IHTt and ([5^ can be recast under the form 



M(Xo) = N\lm^ + ^Xo + Q i^X-''' + ^Xr" ) , (89) 



Xl = Jw? + %Xo \Ax X,y + Br, Xq"" 



(90) 



The sufficient condition to get a closed analytical formula for M{Xq) is to solve analytically This is possible 

if this last equation can be rewritten as a polynomial equation of the fourth degree at most. If A 7^ ?7, it can be 
computed that the following couples will lead to such an analytically solvable equation 

{\^) or (r,,A) = (-^,-1), (0,-1), (1,-1), (-1,-2), (0,-2), {-\,-\) ■ (91) 
Note that the last three cases are only allowed with a nonrelativistic kinematics. 
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B. Equal powers 

1. General case 

The problem actually becomes particularly simple when A = ry. In this case, the following values lead to an 
analytical solution 

A = -1,4 4. 0,1, 2,-2. 4 4 -I 4 4 ,92, 

The last six values are only allowed in the case of a nonrelativistic kinematics. The solution of (|90p reads 



and the mass formula (l89l) becomes 



X,'+' = iA, + B,f(m^ + ^Xo]. (93) 



Mix,) . + + (94) 

A./m2 + %Xo 



Among the values ([M)) . three cases are of obvious physical interest: A = —1, 1, and 2. The case A = 2 corresponds 
to the harmonic oscillator and has been solved in Sec. IIIII The Coulomb problem, i.e. A = — 1, will be specifically 
considered in Sec. IVID while we will discuss the case A = 1 hereafter. Notice that closed mass formulas for any A can 
be obtained in the nonrelativistic and ultrarelativistic limits, as we also show in the following. 



2. Ultrarelativistic and nonrelativistic limits 



In the case m — 0, Eq. allows us to extract the Xq value for arbitrary power 

Q 



Inserting this value in (|94p provides us with the expression of the energy in the ultrarelativistic limit 

A+1 



M{Xo) 



A 



On the contrary, assuming a large mass, mt « M{Xo), it is possible to extract Xq from ([M)) 

Xo = [m{Ax + Bx)]^ . 



(95) 



(96) 



(97) 



Except the Coulomb case (A = —1) which will be treated subsequently, the term is dominant with respect to Xq. 
Expanding the expression of the energy (|94p at lower order in 1 /m leads to the nonrelativistic value of the energy 



MiXo) =mt + 



[Ax + Bx 



(98) 



3. Linear potential 



The case A = 1 corresponds to a linearly rising confining potential. It is of great interest in hadronic physics since 
a linearly rising potential appears to be the best way of modeling the QCD confining interaction within potential 
approaches, see for example Ref. [l^l for more details. For the three-body problem, the one-body potential (term in 
a) corresponds to the so-called Y-junction, while the two-body potential (term in b) corresponds to the so-called A 
approximation. For physical problems, one has the choice of retaining one approximation or the other or both. Here 
for sake of generality we present the general case including both. For a given baryonic situation, one must take A/" = 3. 
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For A = 1, one can write 

Ai+Bi = 

The solution of when A = 1 is given by 



IN , U{U-l) 

-c, with c — a + b\' 



(99) 



(100) 



where F-{Y) is the real root of the equation a;'^ — 3a; — 2Y = 0, that is 



F-iY) 



Y + \/Y^ - 1 



1/3 



Y + \/Y^ - 1 



-1/3 



2 cos I - arccos Y 



for Y > 1, 
for Y < 1, 



and where 



Y 



2Qc 



(101) 



(102) 



As for G_(y), the notation F-{Y) is chosen so that it remains consistent with our previous work Introducing 
this value in the expression of the energy (j94p . a simple calculation allows us to obtain the energy under the form 



M{Xo) = Um 



F-{Y) 
2Y 



F-{Y) + 



F-{Y) 



In the ultrarelativistic limit, where to = 0, one obtains, from (|96p . the very simple relation 

M{Xof = WcQ. 



(103) 



(104) 



Such a linear behavior of the square mass versus the principal quantum number is a well-known fact in hadronic 
physics, where light mesons and baryons are known from experimental data to exhibit Regge trajectories. 
In the nonrelativistic limit, the mass formula (1981) can be recast under the form 



M{Xq) =mt + 



3 fAfQ'^c 



2J2 \ 1/3 



(105) 



The second term is the "binding energy" , which is a positive quantity in case of a positive linear potential. 

In the special case of a linear potential, there is a peculiar relationship between the mass of the two-body problem 
and the mass of the M-hody problem. More precisely, let us call M*^^-' (cr, m, 6, (52) the AFM eigenvalues of the 
Hamiltonian 



i?^^) = (T\Jp^ + rri^ + br 



(106) 



calculated with the natural quantum number Q2 — 2n + 1 + 3/2. Then, using Eq. (|103p and its counterpart for A/" = 2 
in Ref. [§| with the value of the quantum number Q defined in (03]), the mass M'-^\}J' ,m,a,b,Q) of the original 
Hamiltonian is given by 



M(^) (A/-, TO, a, b, Q) = Af (2) (u, m,a + b^J^^^, 



(107) 



VI. BARYONIC-LIKE SYSTEMS 



In this section, we want to discuss a special situation which can have an immediate application for hadronic systems, 
especially baryons (see Sec. IIX|) or glueballs. Let us consider a one-body potential with a linear shape and a two-body 
potential of Coulomb type (A = 1 and ry = — 1 in subsection IV Ap . We noticed that an analytical solution does exist, 
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but one needs to solve a general polynomial of degree 3 and the corresponding solution is rather involved. We prefer 
to give here the solution in the ultrarelativistic limit (m = 0). In this special case the solution of (1901) is quite simple 

Xo = , (108) 



and (|89p can be simplified as follows 



M{Xo) = 2VHy QM - b (^ ^^-^^ ^ . (109) 

Since the argument of the square root must be positive, it is readily checked that, either there exists a maximal 
allowed number of particles for a fixed value of b, cither there exists a maximal allowed value of b for a fixed number 
of particles. The relevance and the accuracy of Eq. (|109l) will be tested in Sec. IIXI 



VII. ATOMIC SYSTEMS 

We turn now our attention to atom-like systems, i.e. systems described by the following Hamiltonian 



+ — a 



^ 1 ^ 1 



(110) 



i<j=i 



For this section, we remain completely general, the only hypothesis being that the parameters a and a should be real 
positive numbers. The obvious application could correspond in a first approximation to M identical electrons, with 
charge e, feeling the attraction of a central static source (the nucleus) and their own repulsion. In this particular 
situation one must make the identification Aira = Me^ and A.'kol — e^. At this stage, we point out that our formalism 
is spin-independent and thus that this model cannot be directly applied to a real physical atom. Moreover, one should 
be very careful in applying the formulas already given. Two cautions are in order: 



• In Eq. (|110p . ri must be the distance between the particle i and the center of mass of the system of Af particles. 
In the atomic interpretation, this center of mass must coincide with the nucleus. This is not always the case 
but this prescription should be valid in case of a spherical atom. 

• Electrons are fermions and the total wave function (space and spin) must be completely antisymmetrical. Since 
the spatial wave function cannot be completely symmetrical, all possible values are not allowed for Q given by 
Eq. ([^ . For closed shell atoms, Qags (SSI) should be used to estimate the ground state energy only. 



Nevertheless, solving Hamiltonian (jllO|) has an intrinsic interest since, to our knowledge, no corresponding analytical 
mass formula is known so far. A general solution can moreover be found with the relativistic kinematics, starting 



from ^ and 



Applying them to the case A = —1, a — a, b — —a and defining 



D 



oM - 



Af [ 2 



1) 



3/2' 



(111) 



the solution of equation ((93)) is given by 



The mass formula ([M)) then reads 



^0 



(112) 



M{Xo) = mtVl -D^- 



(113) 



In the ultrarelativistic limit, this expression gives the value M = 0. This property is well known in the two-particle 
system and prolongates to the A/'-body problem, at least at this level of approximation. Nevertheless, it is quite 
obvious that a system of massless particles cannot be characterized by a (negative) binding energy since the resulting 
mass would be negative. Moreover, with a Coulomb-type potential, the only energy scale of the problem is the particle 
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mass. So, if this mass vanishes, no bound state can exist. In the nonrelativistic regime, one must assume that the D 
quantity is smah in Eq. (jll3p . In this case 

= mf - imtD^. (114) 

It is clear from formula (|113l) that the value of the parameters cannot be arbitrary. One must satisfy the condition 
D < I. Explicitly this means 

aJ^^-a{^ ^%-'^ f\QJ^. (115) 

Let us emphasize that this condition is independent of the mass of the particle [1^. Let us notice that, for Ana — Me^ 
and 47ra = , a maximal value of M exists if Q increases less rapidly that A/"^ . 



VIII. PAIRWISE INTERACTIONS 



In this section, we focus on the case where the Hamiltonian of the system only contains pairwise interactions. This 
is relevant in many physical A/'-body problems. Of course, all the formulae given in the previous sections where we 
considered a mixing of one- and two- body interactions remain valid; it is sufficient to set V{r) = and K{r) = in 
them. But we have additional interesting properties in this case. 



A. Gaussian potential 



We have considered, up to now, potentials which admit an infinite number of bound states. There are however 
many cases in which the number of bound states is finite. One family of such potentials are of the form —ar^ '''' . 
Solutions of the Schrodinger equation with those potentials have been found by using the AFM in Ref. for ij = 1. 
Moreover, recent results have been obtained for semirelativistic iV-boson systems bound by attractive pair potentials 
and in particular for pair exponential potentials [2^. However, Gaussian potentials are sometimes used in A/'-body 
computations, typically in nuclear physics (see for example Ref. [s^l) or Bose-Einstein condensates [3l|. That is why 
we now propose to study the nonrelativistic Hamiltonian 



i—1 — l 



(116) 



Let us perform the change of variables Xi = f3ri, qi — pi/(3. Then we have 



SB. — ^2 Sa 



AT 2 ^ 



i=l 



i<j=l 



ma 
1^' 



(117) 



Only the potential depth g is a relevant parameter of the system. 
The AFM can now be applied to solve i/ga. We have 



In 



(118) 



and consequently the energy spectrum reads, after simplification of Eq. (|76p with — and E — M — mt, 

- 1 . 



In 



(119) 



The equation of minimization of the energy with respect to v is given by 



1/2 , / 



2Q2 



AA(AA- 1)2' 



(120) 
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and its solution reads 

j.o = ge2^«(^) with Y = (121) 

Wq{x) is the principal branch of the Lambert function, defined through Wf){x) e^"^^'' = x (see Ref. for more 
details). After simplification, the energy spectrum is given by 

(Af-l) 4Wo{Y)^ ■ ^ ' 

In physical units, using the rescaling (|117p . the energy spectrum of Hamiltonian (jll6p finally reads 

E,,iY) = ^E,,iY) with Y = - . (123) 

Denoting explicitly by E^'^''{m,a, (3) the previous formula applied to the two-body problem (Af — 2), it is quite easy 
to express the energy £'(^) (m, a, /3) of the A/'-body problem given by Eq. (jl23l) by the duality equation 

An alternative form for the energy of the system is 

^«,„,„,«._^„^.i±p2^ (125, 

This last formula is an upper bound of the exact energy spectrum since g{x) — —ae"^ ^ is such that g" < (see 
Sec. lIVA~3]) . A similar upper bound on the ground state energy has already been obtained in Ref. LS^ by using a trial 
wave function of the form e '^^i<j'^ii . It is interesting to say that, although that trial wave function is similar to 
the harmonic oscillator wave function we obtain in the AFM, both results are inequivalent. The main reason is that, 
in such a variational computation, the scale parameter k is fitted so that the ground state energy is minimal. In the 
AFM, the scale parameter of the wave function is typically proportional to y/rnvQ, where vq is computed from the 
extremum condition (I120p concerning the whole energy spectrum. Thus, even in the ground state case, both results 
do not coincide. 

It is interesting to check that Eq. (I123P has the correct limit /3 — > 0. In the case of a small (3, the potential 
term is approximately equal to X]i<j=i(^<^ + Q^/^^^ij)i ^ harmonic oscillator. As expected, in agreement with 
formula we obtain 



— — — -a + \l^z—Q + <^ir) (126) 

by performing an expansion of formula (|123D in powers of /3. 

It is well-known that exponential-like potentials only admit a finite number of bound states. In particular, the 
energy formula (I123P vanishes for some value of g. These values are called critical coupling constants since they 
correspond to potential strength for which a given bound state appears (or disappears). By solving the equation 
1 -I- 2Wo{Yc) — 0, one finds Yc = —l/{2^/e) and consequently the critical coupling constants read 

For symmetrical spatial wave function, the ground state exists only if g > 3a^|„.^;.^o ~ 9e/2A/'; we obtain in this case 
the relation 



^ (128) 



that has previously been obtained and numerically checked for several exponential- type potentials [32|]. Similarly 

2 

m\n.=U=o = Jf 52|„.=z^=o ' (129) 

indicating that in order to link a A/'-body system one needs a coupling M /2 times smaller than the coupling for a 
two-body problem (s^ . 
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B. Two-body reduction revisited 

In the case of two identical particles, ri ~ r2 = ^'12/2. So, the total interaction is just a function of ri2 and there 
is no need to make a distinction between V and V potentials. For instance, in the case of two identical power-law 
potentials, the total interaction Vt is simply given by 



Vt = sgn(A) + b) 



(130) 



The mass formula coming from the application of the AFM to the A/'-body Hamiltonian ([T]) is defined through 
the equations (|75|) - d7^ . The same procedure can actually be applied to find the spectrum of the two-body-like 
Hamiltonian 



H* = crvV+^ + g V{r), (131) 
where a and g are positive parameters (the * is added to denote the two-body quantities). One finds 

M*{f,*,D*)=<Jfi*o+9V{I{D*o)), (132) 

where the auxiliary fields are solutions of 



A^o = — 
Mo 



2gQ*^D*o 



2.9A*Si^o 



1/4 



1/2 



(133) 
(134) 



Equations (|132p - (|134l) are actually equivalent to (|76| - ((79)) . provided that only pairwise interactions are considered 
and that the following substitutions are made 



9 AA-1 



Mo — Mo, — vq, - — 

(7 Z 



2Q2 



(135) 



Then, if M^'^\cr,g,Q*) is the the mass spectrum of the hamiltonian (|13ip . it follows that the mass spectrum of the 
Hamiltonian 



reads 



or equivalently 
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N -V^'^V N{N -I) 



(136) 



(137) 



(138) 



In these last two expressions, Q is given by (j44p . 



C. Funnel potential 



The funnel potential ar — b/r is one of the most widely used effective QCD potential in hadronic physics. It has 
actually been confirmed by Jattice QCD as being an excellent parameterization of the potential energy between a 
static quark-antiquark pair [27]. But, this potential is also believed to be applicable in light hadron physics, and it is 
of interest to find an analytical mass formula for the Hamiltonian 



Hi 



E 

i<j=i 



ari. 



(139) 
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It has been shown in Ref. @ that the mass spectrum of 

ijf = crVP^T^ + ar- 6/r (140) 

is given by 

Aff ' = 2^Ja((jQ* - b) (141) 
within the AFM in the case m = 0. Applying the duahty relation (jl38p . one is led to the mass formula 

Ml = a^/8JV{JV-l)JVQ- abN^{N - if. (142) 

It is worth mentioning that this last formula can be computed without using a duality relation, but rather by applying 
directly Eqs. ((8T|) and ((82)) . This constitutes a good cross-check of our formulas. 

For completeness, we mention that, in Ref. [lO[, various upper and lower bounds for the mass spectrum of Hamil- 
tonian (|140l) are presented. In our notations, they are mostly of the form 



Mf = min 

r>0 



a\l H J + akir 



(143) 



where P, fci, and are real numbers depending on the approximation scheme which is used (the quantum numbers 
can be contained in P). In the ultrarelativistic limit, the minimization can be performed easily and one gets 



Mf = 2v/afci(crF- fofc_i), 



(144) 



which is formally identical to our formula (jl4ip . This is another consistency check of the AFM since formulas (jl43p 
have been obtained independently by using other methods. 



IX. APPLICATION TO BARYONS 



In order to test the relevance of our method, we apply it for a particular relativistic three-body system: the light 
baryon composed of three massless quarks. In the framework of constituent models, the Hamiltonian for such a system 
is given by [H, IsJ] 

^"=EV^+AEln-i?|-^a5 E (145) 

The dominant interaction is a confinement by three strings with density A meeting at the center of mass. The short- 
range part is given by pairwise interactions of Coulomb nature with a strong coupling constant as- In order to allow 
a more realistic calculation of the correct baryon masses, this Hamiltonian must be completed with terms that can 
be computed in perturbation [SSj. As these contributions represent no interest for the AFM, they are not considered 
in this study. 

Let us note H^^^ the counterpart of once the auxiliary fields are introduced. The eigenstates of iy^-FB 
built with harmonic oscillator states and the eigenvalues are given by the application of Eq. (|109p with Af — 3, a — X 
and b 2a5'/3, 



Mo^™ = ^12A(^B + 3-^j, (146) 

where B is the band number defined as B — 2{ni + 71.2) + {h + h)', see Eq. One can notice the strong degene racy 

due to the particular form of B. This formula generalizes a result obtained for a small value of as in Ref. [33|. 
Following the considerations of Sec. IIV Al These masses are upper bounds of the exact masses. 

In order to test the relevance of this formula, it is necessary to compute accurately eigenvalues of H^. It is possible, 
for instance, to use a variational method relying on the expansion of trial states with a harmonic oscillator basis l35j. 
We can write 

Bmax 

1^-)= E El'^(5''?(^)))' (147) 

B=0 q(B) 



21 



where B characterizes the band number of the basis state and where q summarizes all the quantum numbers of the 
state (which can depend on B). This procedure is specially interesting since the eigenstates of are expanded in 
term of eigenstates of J^afb ^ length scale factor). In practice, a relative accuracy better than 10^** is reached 

with Bniax = 20. Such results are denoted "exact" in the following. 

Before comparing exact masses with formula (jl46p . it is significant to describe the baryon wave functions. Quarks 
are fermions of spin 1/2 with isospin and color degrees of freedom. The global color function is unique and completely 
antisymmetrical. For total spin S = 3/2 (isospin T = 3/2) the spin (isospin) wave function is completely symmetrical. 
The corresponding wave function are of mixed symmetry for S = 1/2 orT = 1/2. Using degenerate eigenstates of 
H^^^, it is always possible to build three-quarks states completely antisymmetrical [2^, [ssj which are characterized 
by the same mass. For instance, when S = T = 3/2 or 5* = T = 1/2, the baryon states possess a spatial wave function 
completely symmetrical. In the real world, the degeneracies are removed by spin- or isospin-dependent operators 
which can be computed as perturbations [3j|. When an eigenstate of is computed for given values of S, T, and 
the total angular momentum L, we determine the band number of its dominant components. It can then be compared 
with the eigenstate of H^^^ with the same spin-isospin quantum numbers and the same band number. 

In Table [H exact masses are compared with the predictions of formula (jl46p for usual values of the parameters 
A and as. The relative error can be large but the quality of the approximation is quite reasonable for a so simple 
formula, and masses are all upper bounds. Due to the strong degeneracy of the harmonic oscillator, all eigenstates 
of if^FB -^j^j^ ^jjg same band number have the same energy. This approximation is better for high values of L. Let 
us note that the relative error can be reduced by a factor around 2 when as — >■ 0. When as = 0, the relative error 
becomes independent of the energy scale factor -\/A. 



TABLE I: Exact eigenmasses of the Hamiltonian (|145|l as a function of the band number B and the total angular momentum 
L for A = 0.2 GeV^ and as = 0.4. The number in brackets is the probability (%) of the component with the band number B 
in the harmonic oscillator expansion ()147|) . These results are compared with the masses given by formulas (1146 1) . H149|) . and 
(|150|) . The number in parenthesis is the relative error (%) with respect to the exact number. All masses are given in GeV. 
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3.584 [86.3] 
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3.625 (1.1) 
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4.253 (14.5) 
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' 0.1) 


3.858 (0.1) 
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3.721 [74.4] 


4.527 (21.7) 
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;io.6) 


3.782 (1.6) 
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3.838 [86.4] 


4.527 (17.9) 
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; 7.2) 


3.895 (1.5) 
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3.966 [93.6] 


4.527 (14.1) 


4.114 


' 3.7) 


4.006 (1.0) 
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4.103 [96.9] 


4.527 (10.3) 


4.114 


' 0.3) 


4.114 (0.3) 



In our previous work about the AFM |4tZll9|, we have shown that it is possible to improve two-body mass formulas 
by changing the structure of the global quantum numbers Q. By fitting another form on exact eigenvalues, a very 
high accuracy can sometimes be reached. In this work, we will proceed differently and will try to use analytical results 
to find the best shape for the mass formula. 

An upper bound of the ground state of can be computed by using the expansion (|147p with just one state for 
which B^i^x = (with harmonic oscillators reduced to Gaussian states). In this special case, a very simple form is 
obtained 



Af^(S,,ax = 0) = — a(3- VSas). (148) 

This result is a better upper bound of the exact ground state than the formula p46p with B = 0. So we can try to 
improve this last formula simply by the changes: 12 — > 32/7r and 2/a/3 — > a/3- This gives 

Mi^FB ^ y^A(i? + 3-V3as). (149) 
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One can see in Table |T] that the masses are greatly improved, but the variational character of the formula (|149p cannot 
longer be guaranteed: some masses are now below the exact ones. Moreover, the problem of the degeneracy remains. 
To cure this situation, we will look at a similar two-body system. 

The masses of mesons composed of two massless quarks can be computed with the two-body equivalent of the 
Hamiltonian (|145p . In this case, for large values of the quantum numbers, it has been shown by the WKB method 
that the mass depend directly of the combination 7rn/2 -|- I [36j . These numbers are in agreement with a result 
obtained numerically in the Sees. 9.1 and 9.3 of Ref. even for small values of n and /. For small quark masses, the 
long-range part of the potential dominates the dynamics of the baryon and we can assume that Eq. (I107P is valid for 
the Hamiltonian (|145p . that is to say that the Q-dependence is the same for two- and A/'-body systems. So, we can 
try to improve formula (jl49D by simply setting 



M2^™ = ^— a(b' + 3- VSas) with B' ^ ^{m + + {h + h), (150) 

which removes the strong degeneracy of the harmonic oscillator behavior. One can see in Table HI that the relative 
error is now around 1%. Despite its simplicity and its nonvariational character, Eq. (|150l) is then a very good mass 
formula for the eigenstates of the Hamiltonian (|145p . It is not sure that the procedure used here to improve the mass 
formula for baryons could work so well for other Hamiltonians. But, this shows that an improvement is possible, at 
least in some particular cases. 



X. CONCLUSIONS 

The auxiliary field method, which is strongly connected with the envelope theory is a powerful tool to compute 
approximate closed analytical solutions of the two-body Schrodinger equation with various types of interactions 0- 
0, Q. The procedure starts with the replacement of an arbitrary potential V{r) by another one V{r) = vP{r) + 
g{v), P{r) being a potential for which analytical eigenenergies can be found, v the auxiliary field and giy) a well- 
defined function of this extra field. With a proper elimination of v, the original potential is recovered. The basic 
idea underlying this method is to consider the auxiliary field as a real number which is eventually eliminated by a 
condition rendering extremal the eigenenergies for the potential V{r). This technique has been recently extended to 
semirelativistic two-body Hamiltonians [§] . 

In this paper, we apply the auxiliary field method to A/'-body problems with nonrelativistic and relativistic kine- 
matics. As this method requires the knowledge of a solvable problem as a starting point, the system of A/'-harmonic 
oscillators is first studied and solved. Some general results are given for A/'-body systems with one-body and two-body 
interactions. If numerical approximate solutions can be easily computed, it is not possible to obtain closed analytical 
formulas for such complicated problems. So, we focus our attention to the case of N identical particles for which 
it is shown that the number of unknown auxiliary fields can be reduced to three or even one. Several problems are 
then studied within this special framework. The case of power-law one-body interactions plus power-law two-body 
interactions is studied and mass formulas are given for atomic-like and baryonic-like systems. Not only the ground 
state is computed but also all excitations. Up to our knowledge, these results for A/'-body problems are obtained 
for the first time. The case of only pairwise interactions deserves a particular attention. In this case, dual formulas 
connecting the A/'-body problem to the 2-body problem with other parameters are obtained. A mass formula for the 
funnel potential is computed, due to its significance in hadronic physics. As an application, a Hamiltonian describing 
a baryon composed of three massless quarks is studied. Accurate numerical eigenvalues are compared with the pre- 
dictions of our method. Results obtained are quite good for the very simple mass formula computed. The formula 
has even been improved by using supplementary analytical results coming from variational and WKB calculations. 
The auxiliary field method yields also analytical approximation of the eigenstates [STj . They are product of harmonic 
oscillator states if the potentials P{x) and P{x) are quadratic ones. Moreover, in the case of identical particles, states 
with a given symmetry can be built by a well defined procedure [26j |. 

Other types of problems could be treated by our method, for instance a chain-like configuration (phonons, gluons 
in a glueball, etc.). For example, the relevance of "polymer chains" of quarks and gluons in a quark-gluon plasma has 
been already discussed in the literature Q. 
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